Neutron Stars as Dark Matter Probes. 
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We examine whether the accretion of dark matter onto neutron stars could ever have any visible 
external effects. Captured dark matter which subsequently annihilates will heat the neutron stars, 
although it seems the effect will be too small to heat close neutron stars at an observable rate whilst 
those at the galactic centre are obscured by dust. Non-annihilating dark matter would accumulate 
at the centre of the neutron star. In a very dense region of dark matter such as that which may 
be found at the centre of the galaxy, a neutron star might accrete enough to cause it to collapse 
within a period of time less than the age of the Universe. We calculate what value of the stable 
dark matter-nucleon cross section would cause this to occur for a large range of masses. 



I. INTRODUCTION 

Observations of the kinematics of self gravitating ob- 
jects such as galaxies and clusters of galaxies consistently 
send us the same message - if we are to believe in Ein- 
stein's theory of gravity on these scales, then there ap- 
pears to be an invisible quantity of dark matter in each of 
these objects which weighs more than the matter we can 
observe. Cosmological observations add weight to this 
hypothesis and tells us that this invisible matter cannot 
consist of baryons, rather it must be a new kind of mat- 
ter which interacts with the rest of the standard model 
rather feebly - dark matter 

The exact coupling and mass of this dark matter is 
not known but has been constrained. One hypothesis is 
that the dark matter annihilates with itself and interacts 
with the rest of the standard model via the weak interac- 
tion. This weakly interacting massive particle (WIMP) 
scenario has gained favour because such particles would 
fall out of equilibrium with the rest of the plasma at such 
a temperature that their relic abundance today would be 
approximately correct to explain the astronomical obser- 
vations. 

Such a scenario also predicts a direct detection signal 
due to the recoil of atoms which are hit by dark matter 
particles, recoils which are being searched for at several 
purpose built experiments (e.g. 0-13 )■ We also expect 
to see signals from the self annihilation of WIMP dark 
matter in regions of the galaxy where the density is large, 
although there are many uncertainties with regards to 
the magnitude of this signal. Neither of these signals has 
yet been detected although international efforts to find 
such signals are intensifying to coincide with the opening 
of the LHC which also may create WIMP dark matter 
particles. 

Since we only understand the thermal history of the 
Universe back to the start of nucleosynthesis, we cannot 
say with any surety whether or not the WIMP scenario 
makes sense. Furthermore there are many other scenarios 
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of dark matter which involve much more massive parti- 
cles or particles which cannot annihilate with themselves 
0, 0]. There is roughly 5-7 times the amount of dark 
matter in the Universe by mass relative to baryonic mat- 
ter. This ratio is rather close to one, a mystery which is 
only solved within the WIMP framework by a happy co- 
incidence. The closeness of these numbers has led some 
researchers to suggest that, like baryons, dark matter also 
possesses a conserved charge and there is an asymmetry 
in this charge in the Universe. If the two asymmetries 
are related then one would require the dark matter mass 
to be approximately 5-7 times the mass of a nucleon. 
This intriguing possibility would be consistent with the 
controversial DAMA experiment Q and the slight hint 
of anomalous noise in the cogent experiment Q. Such a 
dark matter candidate could also have interesting impli- 
cations for solar physics @. 

Since any constraints on the nature of the mass and 
cross section of dark matter particles are interesting, in 
this paper we will consider both of these paradigms and 
see whether or not it is possible to obtain any new con- 
straints from a new angle - namely by considering the 
capture of dark matter by neutron stars. 

The accretion of dark matter onto stellar objects 
has been considered by various groups looking at both 
stars [Tol - [l3 | and compact objects [l5l - [l7j . In partic- 
ular, the ultimate fate of neutron stars which accrete 
non-annihilating dark matter has been discussed before 

[11 El- 

Our aim is to consider the accretion of dark matter 
onto neutron stars in greater detail in order to examine 
whether or not it would ever be possible to either observe 
the heating of a neutron star due to dark matter annihi- 
lation within the object, or the collapse of a neutron star 
which accretes non annihilating dark matter. 

In the next section we will outline our estimate for the 
accretion rate of dark matter onto a neutron star. Then 
we will explain which densities we will be assuming for 
dark matter in the Milky Way. We will then go on to 
work out how hot we can expect a neutron star to get 
simply due to the accretion of dark matter and compare 
this with observations. 

Finally we will look at whether it is at all sensi- 
ble to imagine a situation where the accretion of non- 
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annihilating dark matter onto a neutron star would give 
rise to its subsequent collapse before eoncluding. 



II. CAPTURE RATE OF DARK MATTER 
ONTO NEUTRON STARS 



suggesting the star is unstable. This is of course due to 
the fact that neutron stars are relativistic objects and we 
need to take into account General Relativistic effects to 
calculate the escape velocity properly. 

The space-time geometry inside a static, spherical fluid 
star is 



In this section we will calculate the rate at which DM 
particles will accrete onto neutron stars. The total cap- 
ture rate depends on the density of target nuclei and the 
escape velocity, quantities which both vary throughout 
the star. The exp ression which needs to be calculated is 
the following j20l |: 



C = 47r 



,0 dV 

where the capture rate for a given radius is given by 
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and the various functions within this expression are 
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dy = ^ [erf(6) - erf(a)] , 



^ = m^/Mn, IJL- = (/i - l)/2. 

In the above, is the mass of the DM particle, is 
the ambient DM mass density, A„ is the atomic number 
of the neutron star nuclei, M„ is the nucleus mass, v is 
the DM velocity dispersion, is the star's velocity with 
respect to the zero point of the DM velocity distribution 
and v{r) is the escape velocity at a given radius r inside 
the neutron star (see below); subscript ★ refers to the 
neutron star quantities. 

For neutron stars it is important to note that there 
is a maximum effective cross section - the sum of all the 
cross sections of all the nuclei in the object cannot exceed 
the total surface area of the star since this is obviously 
an absolute upper limit on the total cross section of the 
object. Because of this, cross sections in excess of cr™^^ ~ 
2 X 10""*^ will not increase the probability of capture. We 
will take this into account in what follows. 

The quantities in this equation ^ have various origins 
- V and depend upon the distribution of dark matter 
in the galaxy and we will discuss them later. The cal- 
culation of the escape velocity is more complicated - in 
a normal star we can simply use Newtonian gravity but 
if we were to apply the same simple equations to a neu- 
tron star we would obtain superluminal escape velocities 



ds^ 
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dr^ + r^dfl^, (3) 



where M{r) ~ An p{r')r'^dr' and $ is determined by 
solutions of 
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2 GM(r) 
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where P is the pressure inside the star. Following the 
usual Lagrange method of calculating the escape velocity 
we find that 
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esc 
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V dr ) 



What remains is to obtain the density and pressure as a 
function of radius so that we are able to solve equations 
(121) and (O. This is done using the Oppenheimer-Volkoff 
equations, which are the General relativistic versions of 
the equations of stellar structure. 



dP 

dr 
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We integrate these equations outwards for different cen- 
tral pressures. The initial condition for $(r = 0) is cho- 
sen so that its value at the surface of the star matches its 
solution at large radh $(r > i?Ns) = 5ln(l — 2GM/rc^). 
We then vary the central pressure in order to find the 
maximal mass one can obtain in the given conditions. 
As an equation of state (EOS) for the neutron star 
matter, we consider the unified model developed by 
Pandharipande & Ravenhall [2l| which is based on the 
Friedman-Pandharipande-Skyrme (FPS) EOS. We then 
simply pick a typical solution which gives a profile ap- 
parently shared by the majority of neutron stars [23 |: 
AfNS = 1.44 Mq, i?NS = 10.6 km, pfg"*''^' = 1.4 x 10^® 
kg/m'^. 

Once the structure of the neutron star including the es- 
cape velocity has been obtained in this way we are half 
way to being able to calculate the capture rate by inte- 
grating equation In principle this is not a full calcula- 
tion because while we are calculating the escape velocity 
as a function of radius in a way which makes sense, we 
are not then looking at the effect that curved geodesies 
will have on the capture rate. Nevertheless, we believe 
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that the capture rate calculated with the machinery pre- 
sented here will be accurate within a factor of a few. This 
is an appropriate level of accuracy for this work. 

We also need to make some assumption about the ex- 
pected density of dark matter in the galaxy, which is 
what we shall turn to now. 



III. GALACTIC DENSITY OF DARK MATTER 



profiles since we have not explicitly taken into account 
the effect of baryonic contraction in this work [2^, [2^ . 

We will assume that the velocity dispersion of dark 
matter is 200 km s~^. Sometimes it will be smaller and 
sometimes larger than this depending on the detailed dy- 
namics of the dark matter halo. However, since this is 
a poorly understood subject (30j . we will not attempt to 
model this in any more detail. 
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FIG. 1. Einasto DM density profiles for three different sets 
of parameters: (r_2, a) = (10 kpc, 0.06), (12 kpc, 0.09), (16 
kpc, 0.19) and (20 kpc, 0.53). 



There are various possible methods of obtaining a real- 
istic profile for the density of dark matter in the galaxy. 
First, we choose to use the Einasto profile 



p{r) = p_2e 



(8) 



where p_2 is the DM density at galactic radius r_2 where 
the logarithmic gradient dlnp/dlnr = —2 and a is a pa- 
rameter describing the degree of curvature of the profile. 
We choose this profile because it describes well DM halos 
of various sizes which appear in N-body simulations (23j . 

Several recent studies have estimated the local density 
of dark matter in some detail (23 - [26| . We do not go 
to such lengths and adopt a simpler method, simply en- 
suring that the enclosed mass at the location of the sun 
-^^Mw(^0) = /o^*^ Anr'^ p{r)dr yields the correct Keple- 

rian velocity Wkopi.(i?0) = [GMuwiR)/ Rq^^^ = 220km 
s~^, where Rq is the galactic radius of the Sun. Although 
it neglects baryons, this gives a reasonable density in line 
with other methods and gives us a density of 0.3-0.5 GeV 
cm~'^ of dark matter at the solar radius. One then ob- 
tains a 1-parameter set of solutions in a with the cor- 
responding values of r_2, which lie between 10 and 30 
kpc (sec Fig. [T|). Note that N-body simulations typically 
yield values of a ~ 0.15 — 0.19 for Milky Way size halos 
[27j . We will allow a bit more freedom to adopt steeper 



IV. ANNIHILATING DARK MATTER AND ITS 
CONSEQUENCES 

Now we have the neutron star density and escape ve- 
locity profiles and models for the density of dark matter 
in the Milky Way, we can calculate how much dark mat- 
ter will be captured using equation ([2]). The result as a 
function of radius for the four different density profiles 
can be seen in figure [2] 
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FIG. 2. Dark matter accretion rates vs. galactic radius for 
four different configurations of dark matter (cf. legend in 
figure) and with (r_2, a) = (16 kpc, 0.19). The merging 
of the two first cases is due to the fact that we take into 
account the limiting geometrical cross section of the neutron 
star surface. 



A. Annihilation of dark matter and injection of 

energy 

We next need to analyse what happens to dark mat- 
ter once it has been captured. In order to do this, we 
base our calculations on Kouvaris' work [l6j . adapting 
the formulae to our situation when necessary. For most 
cases of interest where the scattering cross-sections are 
not microscopically small, one can show that the cap- 
tured WIMPs will thermalise relatively quickly, forming a 
roughly Maxwell-Boltzmann distribution in velocity and 
distance around the centre of the neutron star. If the 
particle we consider is an ordinary WIMP which can an- 
nihilate with itself it will do so at a rate determined by 
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the self annihilation cross section and the density of dark 
matter in the star. Following Kouvaris' argument, we as- 
sume that in most interesting situations (where the self 
annihilation cross section corresponds to that required for 
a good relic abundance) , the annihilation rate reaches the 
accretion rate (Fig. [2]) within around 10 million years, 
which is also approximately the time where the DM an- 
nihilation begins to affect the temperature of the neutron 
star. 



B. Cooling of the Neutron Star 

While some energy will be lost to neutrinos, the en- 
ergy generated by the annihilation of most dark matter 
candidates is carried in the main by leptons, quarks and 
photons. We therefore neglect that part lost to neutrinos 
and assume that all the annihilation mass energy goes 
into heating the neutron star. 

Since all the dominant cooling processes involved in 
this situation scale with some positive power of the tem- 
perature except for the WIMP annihilation emissivity, 
the latter will dominate at some point in the life of the 
star. 

The cooling of the neutron star is given by the following 
differential equation: 



dT 
'dt 



EDM 



(9) 



where cy is the heat capacity of the star, edm is the emis- 
sivity (released energy per volume per time) produced by 
the annihilation of the DM when the latter saturates and 
is given by edm = 3Cm^/47ri?^. e^^ is the emissivity due 
to the modified Urea process, which makes the neutron 
star lose energy through neutrino emission by converting 
protons and electrons to neutrons and vice versa, and is 
given by [sij 



(1.2 X 10^ Jm^-'^s 



'^(no) (iO^k) 



(10) 



(where n is the baryon density in the star and uq = 
0.17 fm""^ the baryon density in nuclear matter). Fi- 
nally, accounts for the effective emissivity in photons 
measured in energy over volume and time and is given 

by 



(4/3)^i?3 



1.56 X 10 



13 



I lO^K/ 



Jm ^s 



(11) 

where simply is the rate of heat loss from the surface 
of the neutron star: ~ 'i'n' crT^^^^^^^ (with a the 
Stefan-Boltzmann constant). 



C. Change in outward appearance 
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FIG. 3. Evolution of the internal (above) and surface temper- 
atures of a 1.44 Mq neutron star situated at various galactic 
radii. In the present case, m^= 10 GeV and ao = 1.5 x 10~*^ 

2 
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times. Note that, since when the DM particles annihi- 
late their equilibrium density corresponds to a negligible 
fraction of the total mass, we have neglected the contri- 
bution of the WIMPs to the specific heat. 

In order to determine the thermal evolution of the star 
through its internal and surface temperatures, we solve 
numerically the differential equation ^ . As the temper- 
ature at interesting times (i.e. beyond 10^ years when 
the star starts to cool down) is reasonably insensitive to 
the initial conditions, we arbitrarily set T{t = 0) = 10^ 
K. 

The surface temperature is related to the internal tem- 
perature by the following approximation [32|: 



1/4 



Tsurfacc = (0.87 X 10« K) 
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The DM heats the neutron star at a constant rate 
which will dominate any other thermal process at late 



where gs = GA/nsZ-Rns 
neutron star. 



is the surface gravity on the 
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D. Results for annihilating dark matter 
do neutron stars get? 



how hot 



Changing the position of the neutron star in the galaxy 
changes the accretion rate of dark matter which in turn 
changes the temperature of the neutron star due to in- 
ternal dark matter annihilation. 

Since we arc working with a 1-parametcr set of density 
profiles (see Section UlI)) . our predictions will be model- 
dependent. However, knowing the discrepancy between 
the different DM density distributions gives us a rela- 
tively good picture of the whole space of possibilities and 
since many relations involved in the process are linearly 
related to the density parameter, the reader can scale our 
results up or down as required. 

For every pair of parameters (r-2,a) defining an 
Einasto profile which matches our given constraints, we 
define the properties (todMjCo) of the DM particle we 
want to probe; we then calculate the capture rate for 
our standard neutron star and we finally apply the tem- 
perature gradient given by Eq. ^ in order to find the 
final internal and surface temperatures. In Fig. ([3]), we 
present an example of our results, with = 10 GeV 
and (To = 1.5 x 10"'*^ cm^. 

Using the method described before, we calculate the 
final temperatures of our neutron star at different galac- 
tic radii, starting with tqc = 10~^ pc and going out- 
wards through the halo up to a radius of 50 kpc. For 
every DM density profile and for every type of DM par- 
ticle, we get a specific curve. As we said, however, the 
relation between the various sets of solutions is rather 
straightforward, since the temperature gradient varies as: 
dT/dt cx Cm^. 

In the most favorable cases in terms of DM properties 
and distribution, the highest surface temperatures one 
can possibly obtain lie around 10^ K near the galactic 
centre (Fig. 2]). Given the nature of the neutron stars, 
these final temperatures could result in luminosities of 
the order of 10~^ L^. 



E. Density profiles with central spikes 

In all our calculations so far we have been consider- 
ing straightforward Einasto profiles without taking too 
much care about the complex astrophysical phenomena 
occurring near the galactic centre on a sub-parsec scale. 
The formation of a super massive black hole is though to 
enhance the density in this central region [33| and this 
combined with effects such as self-annihilation, gravita- 
tional scattering of DM particles by stars and capture 
in the supermassive black hole must be included if one 
wishes to obtain a more realistic picture of this central 
region. 

It is interesting to evaluate the possible consequences 
of considering much higher DM densities at small radii. 
To do so, we take as a benchmark the density profiles 
presented in [s^] and we extract a few values in order the 
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FIG. 4. Final surface temperature of the neutron star (at late 
times) for different DM density distributions, with = 10 
GeV and ctq = 1.5 x 10"'*^ cm^ 



determine the evolution of the surface temperature (Fig. 
[5]). One can observe that, in the most extreme situations, 
the neutron star manages to keep its surface temperature 
well above 10^ K at late times if it is situated at less than 
one thousandth of a parsec from the central black hole. 
Since there have been a large number of compact objects 
observed by the Chandra X-ray telescope in this central 
region 35|, this might prove interesting in the future. 
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FIG. 5. Evolution of the surface temperature of a 1.44 Mq 
neutron star situated at various galactic radii (cf. legend in 
figure). The plot is for dark matter with a mass = 10 
GeV and a scattering cross section of ctq — 1.5 x 10~*^. The 
DM densities are deduced from two models: profiles obtained 
by Bertone & Merritt in [sj (BM) and Einasto profiles (E). 



F. Observational Situation 

It is a challenge for astronomers to observe the ther- 
mal emission from neutron stars due to their very small 
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surface area and luminosity. Rather than attempting to 
provide a complete review of this complicated subject we 
will mention a couple of examples. 

Pulsar PSR J0108-1431 is a nearby pulsar located be- 
tween 100 and 200 pc from the Solar System. The spin 
down age suggests that this object is older than 10® yr 
so that its surface will have had a long time to cool. In- 
deed, observations suggest that the surface temperature 
IS a very low Teff < 9 X lO'' K if the distance is 200 pc 
or Teff < 5 X 10^ K if the pulsar is located at a distance 
of 130 pc [3^. Comparison with figure [5] shows that even 
this very low temperature is too high to betray underly- 
ing heating by dark matter, probably by around an order 
of magnitude. 

Pulsar PSR J0437-4715 is a much older object with a 
spin down age of nearly 5 Gyr, yet its surface temperature 
is slightly higher than that of J0108-1431 It seems 
to be challenging to explain this relatively high surface 
temperature (for its age) and it has been s ugg ested that 
the explanation may lie in internal heating [33 • Unfortu- 
nately, it does not seem that this internal heating comes 
from dark matter since in order for this to be the case 
we would require a density of dark matter much larger 
than what is expected in this region of the galaxy (PSR 
J0437-4715 lies only about 130 pc from the Solar Sys- 
tem). Also, if we were to explain the temperature of this 
object using dark matter, we would also expect a similar 
temperature for PSR J0108-1431. 

If it were possible to extend observations in the future 
to lower temperatures and we were able to find a neu- 
tron star with a temperature ^ 10"* K or lower then the 
situation would become much more interesting. 

As we have shown, we can expect higher temperatures 
for neutron stars which lie at the centre of the galaxy 
because of the larger accretion rates in that region. It is 
however difficult to see into the centre of the galaxy, in 
fact the very centre is obscured by dust [13 and as we 
move upwards in frequency through the electromagnetic 
spectrum, we find the centre of the galaxy becomes in- 
visible in the infra-red, only to return into view in the 
x-ray part of the spectrum. Since this rules out the pos- 
sibility of viewing objects with temperatures less than 
many millions of degrees, the observational situation be- 
comes more complicated (not to mention the three orders 
of magnitude drop in luminosity relative to local neutron 
stars simply due to distance). We may be able to see a 
very hot neutron star lying at the galactic centre, but it 
would be necessary to find some reliable way of dating 
the object to a relatively high precision in order to estab- 
lish that this heating is anomalous and caused by dark 
matter. 

From these considerations, it is clear that the conclu- 
sion largely depends on the assumptions made about the 
DM distribution in the Milky Way. Nevertheless, there is 
a small hope that the most extreme configurations might 
be tested through observations. 

Therefore, our chances of using neutron stars as dark 
matter probes through the method described above are 



non-zero but limited. There exists however an alternative 
way of benefiting from the high accretion power of these 
very dense objects if we assume that DM particles do not 
co-annihilate with each other. 



V. NON-ANNIHILATING DARK MATTER 

In the case of non-annihilating DM particles, there 
would be no heating of the neutron star due to anni- 
hilation - the neutron star would simply accrete DM in 
its core. If the amount of dark matter in the core were 
to increase without limit, there would be various pos- 
sible outcomes. In the case of fermionic dark matter 
and in the absence of any pressure due to the exchange 
of gauge bosons, the density would increase until Fermi 
statistics starts to play a role. When this occurs, the 
dark matter core would develop the equation of state of 
non-relativistic degenerate fermionic matter, P cx p^^^. 
In order to calculate the effect of such a degenerate dark 
matter core upon the neutron star, the correct thing to 
do would be to simultaneously solve the Oppenheimer- 
Volkoff equation for the two stars co-located on top of 
each other and therefore both contributing to the gravi- 
tational field relevant for the solution. We have obtained 
such solutions, and have confirmed that when there is a 
large mass contribution due to a degenerate dark mat- 
ter core existing in the star, the normal mass-radius and 
indeed maximum mass of the neutron star that can be 
obtained varies. 

In practise however, the amount of dark matter which 
needs to be added to the core of the neutron star for this 
to happen is typically far in excess of the Chandrasekhar 
mass (Mch) corresponding to the degenerate dark matter 
star for all but the lightest dark matter candidates. At 
this critical point, the degeneracy pressure of dark matter 
would become relativistic and the degenerate dark matter 
core would become unstable. Gravitational collapse of 
the dark matter core would occur, creating a black hole 
at the centre of the neutron star, resulting in the neutron 
matter also becoming unstable. 

The black hole thus created would swallow the neu- 
tron star entirely. This is particularly intriguing possi- 
bility not simply for its inherent drama but also because 
it could in principle provide an additional explanation of 
the unexplained gamma ray bursts observed in the Uni- 
verse other than the coalescence of a neutron star with 
another compact object. It would therefore be interesting 
if enough dark matter could be accreted for reasonable 
values of the relevant parameters - namely the density of 
dark matter in the places where neutron stars reside, the 
mass of the dark matter and its cross section for scatter- 
ing off nuclei. 

Fixing the accretion time available for the neutron star 
to capture DM particles, one can determine, for a given 
m^, the cross-section needed to accrete a mass equiva- 
lent to MchinT-x) which is the amount of accreted dark 
matter required to instigate the collapse of the Neutron 
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star. This procedure allows us to pick out regions on 
the (To — plane where this could occur (Fig. [5]). We 
assume that the Chandrasekhar limit is of the order of 
Mch ~ -^pi/™x consider three accretion times: 

10^, 10^ and 10^° years, corresponding roughly to 0.01%, 
1% and 100% of the age of the Universe respectively. We 
also consider four different densities 10^^, 10^, 10^ and 
0.3 GeV cm-3. The highest density lO" GeV cm'^ is 
the rather extreme limit of the predictions of Bcrtonc 
and Mcrritt for non-annihilating dark matter in a cen- 
tral spike. Since we consider the limiting effective cross 

1e-40 I 1 . 1 . 1 . 1 ' 1 . 1 . 1 



the following way - the typical energy exchanged in a re- 
coil of a dark matter particle with a nucleon goes like 
TUnucv'^ for TO^ 3> TUnuc- Whcu a 10'^ GeV mass dark 
matter particle traveling at 200 km s~^ falls onto a neu- 
tron star its velocity increases to close to the speed of 
light and the energy exchanged in this collision is much 
larger than the total kinetic energy the particle had at 
infinity. Such a particle will lose so much energy in the 
kick it will almost certainly be captured by the neutron 
star [3^. For a much more massive particle this is no 
longer true and it will have a lower probability of being 
captured due to a single scatter, hence the different slope 
at higher energies. 




1e-60 - 



1e-65 I 1 ■ 1 ■ 1 ■ ^ • ^ ■ 1 ■ 1 

100 10000 1e+06 1e+08 1e+10 1e+12 1e+14 
[GeV] 

FIG. 6. Neutron stars immersed for long enough in a high 
enough density of non-annihilating dark matter will eventu- 
ally accumulate the amount of DM particles which correspond 
to the Chandrasekhar limit A/ch and collapse. The red solid 
lines here correspond to the cross section ctq and mass 
required to cause collapse after 10^ years. The green longer 
dashed lines correspond to 10* years and the blue shorter 
dashed lines correspond to 10^" years. The thicker the line 
the higher the density - we plot four lines for each duration of 
time corresponding to (in decreasing thickness and therefore 
density) lO", 10*, lO'' and 0.3 GeV cm'^. Note the least 
dense red line corresponds to such a high dark matter mass it 
lies outside the limits of the plot to the right. The curved line 
corresponds to the latest data from the CDMS experiment to 
place the results in context. The change in the slope of the 
line in the line around 10® GeV is explained in the text. 

section given by the geometrical cross section of the neu- 
tron star {(j^^^ = 2x 10~*^ cm^), there exists a minimum 
value for the DM particle mass below which the fixed ac- 
cretion time is too short for the star to accrete Mch- The 
corresponding limiting (Tq is shown on the graph by the 
upper cut-off at high cross sections. It is possible to see 
that in the most extreme case, namely a neutron star 
immersed in a density of 10^^ GeV cm~^ for 10^° years 
do wc obtain results which are approaching the region of 
interest for direct detection experiments. 

The cross sections required become much smaller for 
higher dark matter masses since much fewer of these par- 
ticles need to be accreted to reach their correspondingly 
smaller Chandrasekhar mass limit. The change in the 
slope of the line at high masses can be understood in 



VI. CONCLUSIONS 

In this work we have investigated whether or not it 
would ever be possible to use the accretion of dark matter 
onto neutron stars in order to understand the properties 
of dark mmatter better. 

We first looked at the effects of annihilating dark mat- 
ter on the temperature of the neutron stars. As can be 
seen in Fig. 21 the highest final surface temperatures 
which could be caused by the heating of neutron stars 
with dark matter lie around 10^ K even in the most opti- 
mistic circumstances. Given the surface area of the neu- 
tron stars, these values would produce luminosities in the 
vicinity of 10^^ Lq with a peak wavelength at about 3 
nm (corresponding frequency ^ 100 PHz). These sources 
would thus radiate mainly in the range between extreme 
ultraviolet (UV) and soft X-rays. Given the important 
absorption due to dust between us and the centre of our 
galaxy and the presence of other luminous X-ray sources 
in this region, we believe that the objects in question 
would prove rather tricky to detect. 

Perhaps more interesting are the constraints on non- 
annihilating dark matter which come from the fact that 
if enough of such dark matter were to accumulate onto 
a neutron star it would form a degenerate star at the 
centre. If this internal star were to get too large, it might 
reach its own Chandrasekhar mass, which is smaller than 
that of the neutron star since the dark matter mass is 
greater than the nucleon mass in most models. In the 
event of the mass of dark matter in the star reaching the 
Chandrasekhar mass of the star, the dark matter would 
lead to the collapse of the neutron star which collected it. 
Such an event might happen for the kind of values of dark 
matter mass and cross section currently being probed 
by direct detection experiments but only in regions of 
extremely high density. On the other hand, for higher 
mass dark matter particles, required cross sections are 
much smaller, since a much smaller mass of dark matter 
would need to be accumulated in order for collapse to 
occur. 

The idea that the accretion of stable dark matter could 
be responsible for the collapse of Neutron stars is very 
exciting, in this paper we have quantified how likely that 
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is. For low mass dark matter particles, it seems extremely 
unlikely. 
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